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(Dated: April 2, 2009) 

It has been argued that a Universe governed by Eddington-Born-Infeld gravity can be compatible 
with current cosmological constraints. The extra fields introduced in this theory can behave both 
as dark matter and dark energy, unifying the dark sector in one coherent framework. We show 
the various roles the extra fields can play in the expansion of the Universe and study the evolu- 
tion of linear perturbations in the various regimes. We find that, as a unified theory of the dark 
sector, Eddington-Born-Infeld gravity will lead to excessive fluctuations in the Cosmic Microwave 
Background on large scales. In the presence of a cosmological constant, however, the extra fields 
can behave as a form of non-particulate dark matter and can lead to a cosmology which is entirely 
compatible with current observations of large scale structure. We discuss the interpretation of this 
form of dark matter and how it can differ from standard, particulate dark matter. 
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I. INTRODUCTION 

There is compelling evidence that baryonic matter in 
the presence of Einstein gravity does not suffice to de- 
scribe the Universe we live in. The most natural and 
popular suggestion is that we are surrounded by a sea 
of massive, non-relativistic particles. Dubbed Cold Dark 
Matter (CDM), it can account for the dynamics of galax- 
ies and clusters and the large scale structure of the Uni- 
verse. Furthermore it can arise in a plethora of extensions 
to the standard model of particle physics and would cur- 
rently exist as a thermal relic from hot era at early times 

Q. 

Alternatives to the CDM scenario have been proposed. 
At the more extreme level, it has been suggested that 
Einstein Gravity is modified, either through higher order 
corrections to the Einstein-Hilbcrt action, or through the 
addition of new gravitational degrees of freedom that af- 
fect the relationship between thegeometric and physical 
nature of the space-time metric 0, Q • The less rad- 
ical proposals typically involve replacing the Cold Dark 
Matter by some non-particulate degree of freedom such 
as a scalar field or a fluid which has an effectively pres- 
sureless equation of state. Within this class of models, 
there have been attempts at resolving both the dark mat- 
ter problem and the dark energyproblem. A notable 
example is that of Chaplygin gas |g, 0] . 

A proposal has been put forward in [8| from a different 
approach. A theory for degenerate metrics is lacking and 
it was suggested that the solution would be the introduc- 
tion of additional dynamical degrees of freedom for the 
space-time connection [f| [l0| . A candidate action, the 
Eddington-Born-Infeld action (EBI) for these degrees of 
freedom was proposed in and it was shown that they 
had unexpected effects: they could mimic the presence 



of dark energy and dark matter in the expansion of the 
Universe and could modify the Newton-Poisson equation, 
leading to flat rotation curves for galaxies. Hence it was 
proposed that the EBI action was a candidate for non- 
particulate dark matter and dark-energy. The Eddington 
action hlj] has also been considered in the context of dark 
energy [ 12|. 

In this paper we wish to study the effect of EBI degrees 
of freedom on the expansion of the Universe and on the 
growth of structure of the Universe. In doing so, we will 
identify the different regimes in the expansion rate and 
how they depend on the parameters in the action and we 
will calculate the effect on the density perturbations and 
the cosmic microwave background (CMB). We will focus 
on two possible uses for the EBI theory, one in which 
the extra degrees of freedom unify the dark sector, as 
proposed in [8( and another in which they co-exist with 
a cosmological constant, playing the role of dark matter. 
As a result we can identify a viable theory of dark matter 
which is competitive with the standard CDM paradigm. 

The paper is structured as follows. In the Section mi 
we display the EBI action and equations of motion and 
rewrite them as a specific case of bigravity or, alterna- 
tively, as a particular bimetric theory, as recently shown 
in [131 ] ; in Section IIII Al we study the dynamics of homo- 
geneous and isotropic solutions to the equations of mo- 
tion (see fl4j for homogeneous but anisotropic solutions) ; 
in Section IIVI we study the growth of linear perturba- 
tions and in Section [V] we calculate the power spectrum 
of density perturbations and anisotropics of the cosmic 
microwave background, allowing us to make a cursory 
comparison with current data; finally in Section I VII we 
discuss our findings. 
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II. THE THEORY: 
EDDINGTON-BORN-INFELD ACTION AS 
BI-GRAVITY OR AS A BI-METRIC THEORY 



The EBI action is 
' ' (fx 



WnG 
+ S m [g] 



-g(R -2A) + -^ v / |g-^K| 



(1) 



where a is a dimcnsionless constant, I a scale, G is New- 
ton's constant, R is the scalar curvature of g^, and S m is 
the matter action. The tensor is the Ricci curvature 
of a connection C" defined in the usual way as 



The connection should not be confused with the 
Christoffcl connection r™„ of the metric g^. Note that 
in the limit in which q — > 0, the action for q^ v reduces to 
the Eddington action [ll[ . The action (fTJ) is a functional 
of g^y and C ^ and is varied with respect to these fields. 

It turns out that there is simpler formulation for this 
theory [l3[ • Define a 2nd cosmological constant 



A- - 



(2) 



and consider the action for two metrics and 



(3) 



As shown in [Dl the action |T]) is fully equivalent to ^ . 

0% is r 
usual metricity relation, 



The connection Clg is related to the metric q^ v by the 



C£, = 5(«r'i 



a/3 



+ dyq^fj - dpq^) . 



where (q is the inverse of q^„ such that 

(<r i r"<3w = ^ v 

and K ee (q-^K^. 

Theories of bigravity have been proposed in a number 
of contexts: as spin-2 theory of the strong interaction 
EH ! as a full non- linear extension of the Fierz-Pauli 
theory of massive gravity, and more rece ntly as an effec- 
tive theory of interacting brane- worlds [17] . A number 
of examples of bigravity theories have been studied in 
detail [HI, [Til, H3, EH , in particular in terms of their con- 
sistency, asymptotic behavior and the global dynamics 
of isotropic and homogeneous space times. It should be 
remarked that in the context of bigravity theories, EBI 
turns out to be the simplest theory with a minimal inter- 
action between both sectors. In this paper we study in 
detail the cosmological dynamics both of the background 



and at the perturbative level and hence extract useful 
hints of what one might expect from more general classes 
of bigravity theories. A class of bi-measure theories have 
been considered in [22| and references therein. 

There is yet another point of view one can take of this 
theory. If one looks at the action as it is presented in 
equation^ i.e. a theory of two metrics, one of them, g^, 
quite clearly couples to the matter fields and has physical 
significance- it is this metric that defines how clocks and 
rulers respond- and hence we can call it a "physical" met- 
ric. The other metric, q^ satisfies the Einstein-Hilbert 
action and couples to the rest of the world through it's 
interaction with the physical metric. If we interpret q^ u 
to be the metric of space time- we can dub it the "geo- 
metric" metric- we then have a bona-fide bimetric theory 
of gravity. This is entirely akin to the approach in the 
Tensor- Scalar- Vector theory of gravity [4| and can give 
us an intriguing interpretation of roles of the different 
fields. 

The field equations which are found from cither the 
original EBI or the bigravity (or bimetric) action are the 
Einstein equations for g^ v 

G\ = BnGT»„ - A*"„ - LJ1 {q- l f a g av (4) 

and the Einstein equations for q^ v 



(5) 



where Q M „ = K^ v — ^K5^ v is the Einstein tensor of 
q^u . These are the complete set of equations with which 
we can study the dynamics of the EBI action. (Note 
that tracing (O a simplified equation is obtained K^ v = 



III. COSMOLOGICAL DYNAMICS 

A. FLRW equations 

We now focus on the dynamics of homogeneous and 
isotropic metrics in EBI gravity and will restrict ourselves 
to spatially flat metrics so that the line element is given 
by 



ds 2 



-dt 2 



where t is physical time, x are spatial coordinates and 7y 
is the metric of a flat hypersurface. Note that from the 
bimetric point of view, this makes sense- all observables 
will depend on the physical metric, g^ v . The vanishing of 
the Lie derivative for all Killing vectors of the spacetime 
gives q^ v such that, 



9oo 



-X" 



Hi 



(6) 



The functions X, Y parameterize the metrics compatible 
with the background symmetries. 
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The Friedmann equation for this cosmology is 

3H 2 = 8nG(p E + p f ) (7) 
where H — -, the EBI density is given by 



- Y3 
pE ~ 8wG£ 2 Xa 3 



(8) 



which is always positive, and pf is the energy density 
in all the remaining fluids (including the cosmological 
constant, A). The Raychaudhuri equation becomes 



H 2 = 8ttG(P e + P 



f) 



where Pe = — XY/8nG£ 2 a and Pf is the pressure in 
all the remaining fluids (including the cosmological con- 
stant, A). The remaining held equations are then 



3t 



Y 2 
6£ 2 — 
ut y 2 

Y Y X 

Y ~ YX 



Y 2 



1 + aX 2 



(9) 



The EBI degrees of freedom behave as a fluid. We can 
trade the variables X and Y in terms of the fluid density 
pE and equation of state parameter we given by 



a 2 X 2 
Y 2 



w E 



(10) 



Eliminating the coordinate time t using the Friedman 
equation (|7|), we find that the new variables evolve as a 
function of ln(a) as 



p' E = -3(! + w e )pe 



(11) 



and 
w' E = 2w E 



1 + 3we 



+ \U(-w E ) 3 / 2 n E a-2 



(1 + 3w E )pe 



Pc 



(12) 



where the relative density of a species " i" (including EBI) 
is as usual Cli = Bi. for the critical density p c = pE + pf, 
and where we have defined pg = [87rG^ 2 ] -1 . Clearly the 
above equation is inconsistent for we > 0; in fact we is 
bounded from above by the condition 



16TrG£ 2 (-w E f /2 pEa - 3w E > 1. 



B. Analytical approximations 



(13) 



The system ([7]), (fTT|) and (|12|) can be solved numeri- 
cally but we can extract some analytical results by ex- 
amining its approximate behaviour. 



The first clear case that one can see, is that when 
we ~ 0, the equation of state we will evolve slowly 
and the EBI field will behave as cold dark matter. This 
is independent of whether the EBI field is dominating 
the background dynamics or not, and is therefore valid 
throughout the entire history of the universe (e.g. dur- 
ing radiation, matter and possible cosmological constant 
eras) provided \we\ is small. In this case we ~ wqo 2 
where wq < is an initial condition. Since We is propor- 
tional to a 2 , the cold dark matter behaviour is unstable, 
and is bound to end when we is sufficiently driven away 
from zero and becomes 0(1). At this point, the sub- 
sequent behaviour of the EBI field, depends on various 
factors which we analyse below on a case by case basis. 



1. Case : ps ^> Pf 

Let us first consider the case where pe 3> Pf i.e the 
EBI field is driving the background dynamics. Apart 
from the dark matter behaviour (we ~ 0) which as we 
have discussed above can always be realized, we uncover 
two more phases. The first is a constant-w phase such 
that we — w c which solves the equation 



(1 + 3w c ) 4 + 16a 2 < = 



(14) 



This constant w depends only on the parameter a. For 
a = we have that w c = ~ \ while for a = 1 we have 
that w c = — 1. This means that — 1 < w c < — i for 
< a < 1 and w c < — 1 for a > 1. Indeed in this case 
phantom behaviour is possible by allowing a > 1. 

When a < 1 (i.e. —1 < w c < — |), the above constant- 
ly phase is unstable, simply because pe eventually drops 
and approaches p£. If this happens before the fluid pf 
becomes dominant (which could happen if pf is a cos- 
mological constant), then the constant-w phase ends and 
the EBI fluid now behaves like a cosmological constant, 
i.e. we = — 1. Eq. (fT2|) . then gives that this cosmological 
constant is given by 



PE 



1 



a 



(15) 



It turns out that this deSitter phase is stable under ho- 
mogeneous time-dependent perturbations (although not 
under inhomogeneous perturbations; see relevant section 
below). We postpone the stability analysis until case 3. 



2. Case : pe = Ppf and pf 7^ pA (tracking phase) 

An interesting case emerges if pE is neither negligible 
nor dominant but rather is assumed to track the fluid. 
This is possible provided pE — Ppf where (3 is a propor- 
tionality constant related to the EBI relative density as 
VL E = (3 /(l+ (5). Eq HU) then tells us that this is possible 
iff w < — i, in which case we get 



(l + 3w) 2 



A(-wf/ 2 a- (\ + 3wf 



4 



which gives fl E = 4 \J w) % a . 

Since there are no interesting fluids with equation of 
state w < —1/3 (apart from a cosmological constant, 
treated below) , this case is not of much relevance. 

3. Case : pE = Ppa (cosmological constant tracking phase) 



Therefore the approach to deSitter is critically damped 
for a — a c = g^Ej^ while it is underdamped for 
a < a c and overdamped for a > a c . Furthermore, we 
have that -f < 1 - |(1 - a) Pl ~ PK < 1 and so the 

3 3V I pi - apK 

underdamped solutions are always decaying. Finally it is 
possible to have only the overdamped solutions by choos- 
ing p A < p e < §p A . 



In the limit where w — —1 we find (3 = l/(a — 1) and 
therefore pe — Pa/ {a — 1), which is valid only for a > 1. 

We can explore this limit further; if pe is constant and 
pi non-negligible then we have that the EBI density is 
given by 



Pe 



Pi - PA 
1 — a 



(16) 



while the effective cosmological constant such that 3H 2 



8nGp 



(off) 



IS 



4. Case : p E < pf / pA 

We now pass to the regime where pe -C pj, i.e. the 
EBI field is subdominant, and the cosmological dynamics 
are driven by some fluid pf which is not a cosmological 
constant. Consistency requires that pi <C pf (otherwise 
this case cannot be realized), and we therefore get that 



w' E w 2w E {l + 3we) 



(18) 



(eff) Pe - apA 
p\ ' = Pe + Pa = —rz 



(17) 



(off) 

Pa 



In the limit in which £ — ► 00 we recover the previous case. 
One should further impose the conditions pe > and 
> 0. For < a < 1, a necessary and sufficient 
condition for this to hold, is that pi > pa (regardless 
of the sign of pa), while for a > 1 we need < pi < 
Pa- This second subcase cannot be realized (see case-5 
below). For a < 1, taking the limit pa — > takes us back 
to case-1. 

The negative sign appearing in the expressions above 
is quite misleading, and one could think that it might 
be possible to cancel the effective cosmological constant 
to sufficiently small values. This is clearly impossible 
for pa > simply because by virtue of ([H]) wc also 
have pe > 0. It is also impossible for pa < since 
again because of we need a < 1 which implies that 

Pa ^ ^ P E > Ipa! • Thus we cannot have cancellation of 
the cosmological constant. 

We now perform stability analysis (as mentioned in 
case-1) of this deSitter phase (for which a < 1). Let 
PE = T-rfMl + d) and w E = -1 + e 2 with e x > and 
€2 > 0. Perturbing (fTTj) and (fT2]) to linear order we find 



' = -3e 2 



Hence we find that if the EBI-field is subdominant, the 
above equations lead to two possible behaviours : the EBI 
fluid behaves either as cold dark matter ifO<W£<C— 1, 
i.e. very close to zero, or as curvature if we ~ — i. 



5. Case : pE *C pf = pa 

The final case we consider is when pe is negligible but 
now pf = pa, i.e. the background fluid which drives the 
dynamics is the bare cosmological constant A. Here we 
find a new regime such that 



w E = we 



3 PA 



(19) 



in addition to the cold dark matter regime which can 
still be realized. Note that we recover the curvature like 
behaviour if pi -C pa while we = — 1 if Pa = Pe — 
(8WG1 2 )- 1 . 

One may wonder whether phantom behaviour such 
that we < — 1 can be realized in this case, by choosing 
I such that 8itG£ 2 pa < 1. This turns out to be impossi- 
ble: when pa becomes smaller than the threshold value 
8itG£ 2 pa = 1, this takes us back to case 3. Decreasing I 
further (or decreasing pa) eventually leads to case 1. 



and 



/ on i Pe - Pa 
e 2 =2(1 ~ a) ei - 3e 2 



Pe - oipA 



which combine to give 



" i o ' i c/i \ Pe~ Pa n 
ei + dei +6(1 — a) 61 = 

Pe - apA 

The normal modes are e" ln a where 



n = — 
2 



1 /-, 8 p e - p A 

3 pi — apA 



C. Realistic model building for the background 
dynamics 

Having analyzed the different possible behaviours of 
the EBI field in various cases above, we now turn to re- 
alistic model building. 

We will use units in Mpc which is the standard in pop- 



ular Boltzmann solvers such as CMBfast [23j, CAMB [24 [ 
and CMBeasy [25] . In these units we have that the Hub- 
ble constant today is Hq = 3.34 x 10~ 4 h Mpc' 1 , with 
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FIG. 1: Upper panel : The locus of fixed to and angular 
diameter distance to recombination, by varying ujl and wq 
keeping all other parameters fixed. The value uja = 0.36 
corresponds to the WMAP5 best fit model. The other EBI 
parameters are 1 — a = 10~ 6 and £ = fO 9 Mpc. Lower panel 
: A similar locus curve, only now we vary a and wo keeping 
all other parameters fixed and in particular uja = 0. The 
alphabetical labels are explained and discussed in the main 
text. 



h ~ 0.6 — 0.8. For numerics we can absorb 8irG into the 
definition of densities. The total fluid density is thus 



SirGpf = 3.34 x 10" 



UJ b 

-^T +LU A 

a 6 



Mpc 



(20) 



where uj r = 4.16 x 10~ 5 (for CMB temperature of 2.726iir 
and three species of massless neutrinos), uib ~ 0.018 — 
0.023 given by nucleosynthesis, and uja ~ — 0.5. 

Turning to the EBI field, we need to set its initial den- 
sity f>E,in and initial equation of state parameter WE.in 
at the initial scale factor aj. To do this we require 
that the initial condition for pE is such that it would 
give rise to an equivalent CDM density in the past. In 
other words requiring that the equivalent CDM density 
today would be 8TrGp c = 3.34 x 10~ 7 uj c Mpc~ 2 , with 
uj c ~ 0.11, we extrapolate this to the initial scale fac- 
tor a,- and set the initial condition for the EBI den- 



sity as 8nGp E in = 3.34 x 10 



-7 



x a„ 



uje Mpc 2 , with 



w E ~ 0.08 -0.13. 

For setting the initial condition for w E , we require that 
the EBI field behaves as CDM all the way up to at least 
a ~ 0.1 where we starts to become O(l). Since in the 
CDM phase we ~ —w a 2 , we set the initial condition for 
We as WE,in = —wocq, by specifying a positive parameter 

Wq. 

The background model is thus completely determined 
by six parameters : the initial conditions UJb, u>e and 
wo as well as the two parameters £ and a (for the fixed 



radiation density discussed above; massive neutrinos can 
easily be accommodated in the usual way but we refrain 
to discuss it here for reasons of simplicity). On top of 
specifying these parameters one has to make sure that 
the inequality constraint (fl~3|) is obeyed. In the light of 
setting up initial conditions as we have just described the 
inequality becomes 



a 



( 1 



\Mpc 



3/2 



ue > 1-5 x 10 6 



(21) 



1. The AEBI model 

The simplest possibility is when the EBI field is chosen 
to act as CDM all the way, and is not responsible for the 
accelerated expansion of the universe which is due to the 
bare cosmological constant A. In this case, in order for 
the EBI field not to deviate from the CDM track, we 
must have that wo < 10~ 4 . Moreover the parameters I 
and a do not have any role in the background. Thus we 
are down to three parameters : u>b, uja and uje, the same 
as in the standard ACDM model. We call this model the 
AEBI model. 

Choosing UJb — 0.023, uja — 0.36 and uje = 0.114 we 
have a background evolution identical to the WMAP5 
best fit ACDM model. It is therefore not distinguish- 
able from the ACDM model using, for example, type-la 
supernovae data [26| . 



2. The general EBI model 

By lowering the cosmological constant wa gradually 
to zero, we should start compensating by having the EBI 
field to play a role in the accelerated expansion of the 
universe. 

Lowering A changes the angular diameter distance to 
the surface of last scatter and hence shifts the position 
of the acoustic peaks in the CMB. Since the peaks are 
very tightly constrained, we must change one further pa- 
rameter to make up for the change induced by varying 
A. As a 2nd parameter (for illustration) we choose to 
vary wq. The upper panel of figure [1] shows the locus 
curve for which both the angular diameter distance and 
conformal time today To, remain constant, in the wq and 
uja plane. The point uja = 0.36 and wq = 6 x 10~ 8 . 
We call this model " Model A" . The parameters a and 
£ are chosen to be |1 — a\ = 10 -6 and £ = 10 9 Mpc for 
this model. Increasing £ or lowering wq still gives ac- 
ceptable AEBI models as the dependence on lowering wq 
or increasing £ is very week and for most cases does not 
produce any observable result. As we lower uja towards 
zero, we keep increasing wq, along the displayed curve to 
make sure that the angular diameter distance to recombi- 
nation stays constant. When uja — exactly, wq « 0.845 
(this depends on our choice of uib and uje)- We call this 
"Model C". For clarity we also consider a "Model B" 
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FIG. 2: Four models from the locus curve of Fig[T] (ua,»o) 
= (0.36,6 x 10" 8 ) (upper left), (0.2,0.43) (upper right), 
(0.1,0.65) (lower left), (0.,0.84) (lower right). The curves 
are : EBI (solid), A (dashed), baryons (dotted) and radiation 
(long-dashed). Note that in the last case (lower right), the 
final state of the EBI field is an approximate cosmological 
constant. The actual phase is the constant-w phase with w 
being extremely close to —1 due to 1 — a = 10 -6 . 



1000 
100 
10 

1 

0.1 

>" 0.01 
0.001 
x 0.0001 

io- 5 
io- e 
io-' 
io- 8 
io- B 
io-'° 
io-" 
io-' 2 

: o» 
:o> 

10" 

000 
100 
1 

1 

0.1 

>< 0.01 
- 0.001 
X 0.0001 

io- 5 
io-« 
io- 7 
io-» 
io- 9 
io-'° 

10 



10 



o.oooa.ooio.oi o. 



10 0.0003.0010.01 0.1 1 10 



FIG. 3: Same four models as Fig[2] They are from the lo- 
cus curve of Fig[T] (w A , w ) = (0.36,6 x IO -8 ) (upper left), 
(0.2,0.43) (upper right), (0.1,0.65) (lower left), (0.,0.84) 
(lower right). The curves are : X (solid), Y (dashed), 

■^■(dotted) and y (long-dashed). 



for which u>a — 0.2 and wq « 0.43. Figure [2] shows 
the evolution of the EBI energy density (solid), radiation 
density (long dash), baryons (dotted) and cosmological 
constant (dash) for four models along this curve. The up- 
per left panel is Model A, the upper right panel is model 
B with wa = 0.2 and wq ~ 0.43, the lower right panel is a 
model with lu\ — 0.1 and wq ~ 0.65, and finally the lower 
right panel is model C. Figure [3] exhibits the variables X 
(solid), Y (dashed), ^ (dotted) and y (long-dashed) for 
the same set of models. Observe that during the time for 
which the EBI field is like CDM, both X and Y are ap- 
proximately constant (very slowly varying) , while during 
the we ~ —1 phase, X is still an approximate constant 
while Y is varying. 

Once we reach the point u\ = (model C) we can 
start investigating the effect of changing a. Once again 
we keep changing wq in order to compensate and keep the 
angular diameter distance the same. Along this line we 
consider two further models : Model D has a = 0.3 and 
wo w 1.166, while Model E has a = 2 and wq w 0.632. 
Remember that a is connected with the constant we 
phace. In particular Model D has we ~ —0.493 during 
the acceleration era while Model E gives rise to phantom 
behaviour with we ~ —1.795. 

In principle we can further investigate varying I. We 
find however that changing I does not lead to any inter- 
esting new behaviour. If we compensate the variation of 
£ with wq, then wq must be increased, and this leads to 
a similar effect as decreasing a. 

Let us also note that the general EBI family of models 
would give a background evolution that deviates from the 
ACDM model. They can therefore be distinguished from 
ACDM using, for example, type-la supernovae data [Hj]. 
However, as we discuss further below, we do not find 
any parameter space allowed (apart from the AEBI spe- 
cial case) when we consider the Cosmic Microwave Back- 
ground angular power spectrum observations. It is there- 
fore of little significance to try to constrain such models 
with the supernovae data. 



D. Summary of the background evolution 

In this section we have mapped out various possibilities 
for the background evolution. As claimed in [9], it is 
possible to construct a theory in which the EBI field plays 
the dual role of both dark matter and dark energy. We 
have generalized the results in Q by uncovering a third 
phase of the EBI field where it has a constant-w equation 
of state which interpolates between the CDM phase and 
the cosmological constant phase. Thus the EBI field is 
a unified model very similar to to Chaplygin gas. We 
stress however that contrary to the Chaplygin gas, we is 
an independent dynamical degree of freedom, and thus 
the equation of state of EBI is not rigid. This leads to 
even richer dynamics in the perturbations as we show in 
the next section. 

We have specified the requirements on the initial con- 
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ditions of the EBI field as well as the two parameters £ 
and a, in order to have a background evolution that is 
compatible with the standard paradigm. This gave us 
the simple AEBI model, for which the EBI field replaces 
CDM but is not responsible for the accelerated expan- 
sion which is still due to the cosmological constant A. 
By gradually lowering A to zero, and compensating with 
the EBI field, one can have the initial condition wq as 
well as the parameters £ and a to play a role, leading 
to the mixed EBI model, where the effective cosmologi- 
cal constant recieves a contribution from A and the EBI 
field. When A = we get the plain EBI model. 



equations are 



2k 2 ® - 6-$ - <J 



2$ + 2-* 

a 



6$ + 6- 

a 



2$ + # 



. a a 



2k 2 (<S> 

2 ' 



* = 



a 



87rGa 2 ^ ptfi (22) 

i 

8TrGa 2 ^2pi&i (23) 

i 

*) 

247rGa 2 p^ (24) 

i 

S-KGa 2 ^2piS t (25) 



IV. EVOLUTION OF INHOMOGENEITIES 



As we have seen in previous sections, the EBI back- 
ground field can behave as pressureless matter and as 
cosmological constant, and there are transitions between 
these two phases. 

After displaying the equations of motion for linearized 
fluctuations, in this section we show via analytical ap- 
proximations that the matter phase is consistent with 
current observations. On the contrary, fluctuations on 
the acceleration phase show an unacceptable growth. Fi- 
nally, we show an analytical series in l/£ providing a 
systematic way to isolate the matter phase, as the zero 
order approximation, leaving an evolution indistinguish- 
able from ACDM. Corrections in powers l/£ can then be 
computed to any desired order. These will be reported 
elsewhere. 



where the index "i", runs over all fluids, including the 
EBI field. Here S is the fluid density fluctuation, O is the 
fluid momentum divergence, n is the pressure fluctuation 
and S is the fluid shear. We can combine them to find a 
Newton-Poisson-like equation of the form: 



2k 2 ® = 8nGa 2 Y^ 



PiSi + 3-piQi 



(26) 



The EBI density contrast Se, velocity perturbation 
&e, relative pressure perturbation He and shear per- 
turbation, Se, are defined as linear combinations of the 
EBI metric variables S, (3, p and \ : 



d E = *-S + 3($- X ) 
Oe = -P 

U E = w B (S-* + $-x) 
S = —weP 



(27) 
(28) 
(29) 
(30) 



The evolution equations for the above fluid variables 
are found to be 



A. The equations of motion 

For the purpose of studying large scale structure, i.e. 
density perturbations, we will focus on scalar pertur- 
bations in the conformal Newtonian gauge in confor- 
mal time, r, such that goo — — a 2 (I + 2^f), got = 0, 
gij = a 2 (I — 2$)7y . We have that Vi is the covariant 
derivative on the hypersurface such that Vijjk = 0, and 
define D i3 — V^V ' j — ^7^ V 2 . 

The tensor field q^ is perturbed as qoo = —a 2 X 2 {l + 
2H), q w = -Y 2 V z /3, qij = Y 2 [(1 - 2 X Hj + D l3 p}. No- 
tice that the g-metric has four scalar modes, namely S, 
j3, x an d v, as there is no gauge freedom left to set any 
of them to zero. We also find it convenient to define 

17 _ dlnY 
dr 

As in the homogeneous case, it turns out that the EBI 
field can be cast as a generalized fluid, in the framework 
of a generalized dark matter model [27|. The Einstein 



6 E = -k 2 E + 3(l + w E )<i> + 3-(wE6E~TlE) (31) 

a 

e E = -(3wE-l)6E + (l + WE)y-lk 2 SE + n E (32) 
a 6 



, a weM 2 

4Z + 2(1 + 3 Wb ) f— 

a iZ 



2we 

w e 
~2Z 
1 

~2Z 



1 



w E 



we 



3a 2 



9-2^$ 
Z 



2£ 2 k 2 

a 2 {3w E -l) 
2£ 2 k 2 

3a 2 (l + w E ) 
f 2£ 2 k 2 



Se 
He 



cl 2 (1 + we) k 2 , h 



-WE 

l l- 2 

-g« w E 



a {3w E -l)a 2 , k 2 

- Z ~ i a WE+ U 2 Z + 3Z WE 



+w E 



9k 2 
4fc 2 



4Z* - — w E <f> + (1 - 3w E )® 
0Z1 



(33) 
U E 
Se 

(34) 
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The remaining fluids can be described in the usual way 
using conservation of energy and momentum. 



B. EBI as CDM 

In this section we show how the EBI field can behave 
as CDM even at the perturbative level. First notice that 
after setting He and Se to zero in all of the fluctuation 
equations, the remaining two variables, namely 5e and 
e will obey equations which are the same as for a CDM 
fluid, provided the background equation of state param- 
eter we is also very close to zero and the EBI field is in 
the background CDM phase. 

We now show that if we ~ and the pressure pertur- 
bation He and shear Se are initially chosen to be zero, 
then they will remain arbitrarily small, and the EBI field 
will behave as CDM even at the fluctuation level. 

During this phase we get that Z as 
aH^aflEi—WE) 3 ^ 4 , hence we can set we — * 0, 
Z — ► and we/Z — > 0. The IIe equation then becomes 



dint 



n, 



(35) 



Since it is not sourced by any other variable in this limit, 
if we set He = initially, it will stay zero. Thus imposing 
He = in the shear equation and taking the same limit 
above we get 



dS E 
dint 



2S r 



(36) 



which has solution Se = Sqo 2 . 

Thus for very small initial relative pressure perturba- 
tion and shear, i.e. n ~ and Se — 0, the EBI field will 
behave exactly as cold dark matter at the perturbative 
level, i.e. the equations for 5e and Qe would read 



-fc 2 6i; + 3<i> 



®e = — + * 
a 



(37) 
(38) 



respectively. 

This means that if \we\ <C 1 throughout the entire 
history of the universe until today (such as the AEBI 
model-A), we would expect any observable to be com- 
pletely indistinguishable between an EBI dark matter 
model and a standard dark matter model. Note also that 
we have assumed that IIo ~ and Se — to obtain exact 
CDM-like behaviour. But if we were to include a small 
amount of pressure perturbation and shear in the initial 
conditions on very small scales, this might lead to differ- 
ences with the CDM which might be observable at the 
cores of galaxies, clusters and in the small scale structure 
of gravitating bodies. Indeed, differences vis-a-vis CDM 
are expected to occur once a system enters the non-linear 
regime where EEq and Se — are inevitably sourced. 



C. EBI acceleration era 

The evolution of perturbations during a regime of ac- 
celerated expansion is more intriguing. Assume that 
A = and neglect baryons and other components. As 
shown in 0] the equations of motion for EBI gravity have 
an exact de-Sitter solution, 



o(r) 



TO 



X(r) = 



1 



Y(r) = 



TO - T 



This field solves equations and © • tq is an arbitrary 
integration constant fixed by the initial conditions, a — > 

00 aS T — ► Tq . 

To check whether this solution can represent or not the 
accelerated era of the Universe we study fluctuations on 
this background. This is an straightforward exercise and 
we summarize here the main results. Using the equations 
of motion all EBI functions x(r), /3(t), H(t), /j,(t) can be 
written as functions of the Newton potentials 3>(t) and 
^(r) by algebraic expressions. We are left we two coupled 
second order differential equations for $(t) and ^(r). 
Interestingly, the combination $(r)+^ r (r) decouples and 
satisfy a Bessel-like equation with the exact solution, 



$(r) + *(t) = 



Q J u {k(T - t)) + b Y„(k(T - r)) 



where 



v = -\/24a - 39. 
2 



(39) 



For any value of a this function diverges as r — > To mak- 
ing de Sitter space unstable. 

We would like to stress that this conclusion may change 
when introducing extra "Fierz-Pauli" couplings to the 
EBI action. See UM for a recent discussion. 



D. l/£ and the AEBI theory 

In this paragraph we develop an analytical method to 
isolate the matter phase for the EBI field. Numerical 
analysis show that this occur for £ and Y large. We then 
attempt to use l/£, the coupling between both metrics, 
as a perturbative parameter. 

Consider the following Frobenious type series for the 
background functions a(t), X(t), Y(t), 

a(r) = a (r) + ^(r) + ^a 2 (r) + ■ ■ • 
X(r) = Xo(t) + ±X 1 (t) + ±X 2 (t) + --- 

Y(r) = ^/^y o(r ) + Ir 1 ( T ) + ^y 2(T ) + ...^(40) 

We have included a positive power of £ in Y for two rea- 
sons. First, recall the EBI background density has the 



9 



form pe ~ ~p~Ta5- Thus, if Y scales as Y ~ £ 2 / 3 , then at 
order zero in 1/2 there will be a finite contribution to the 
Friedmann equation from the EBI field, which turns out 
to be dark matter. The interesting observation is that 
this series also provides the right equations for all other 
variables including fluctuations. 

There is another reason to include a growing £ term in 
the background. If the metric functions do not depend 
on £, then in the limit of large I the metrics g^ v and 
q^ v become decoupled. This is clearly seen from the bi- 
gravity action The decoupled system describes two 
massless gravitons and is a theory with a different num- 
ber of degrees of freedom. Switching on the interaction 
term, proportional to 1 /£, is in this sense a discontinuous 
change to the theory. On the contrary, for a series of the 
form (|40|) the two metrics never decouple. 

We first discuss the background equations and then 
the fluctuations. Plugging (|4"0)) in the background equa- 
tions we find consistent equations for the coefficients 
ci,-(t), Xj(r), Yi(r) order by order in I. 

At order zero one finds that Xq and lo must be con- 
stants while ao( r ) satisfies the Friedman equation 



^o 3 



1 



1 



4(t) 3Xo4(t) ' 3' 



-A 



(41) 



Thus, by choosing Yq /Xq and A appropriately we find 
a background evolution which is, at this order, indistin- 
guishable from ACDM. 

Corrections 1/ £ can be computed order by order. Since 
£ is so far arbitrary we can make it as large as necessary 
in order to suppress l/£ corrections. For completeness 
we display the first order equations 



dY t 



Y 



6 (^L) -2aX 2 Y 2 a (r)-Y 2 a (r) 2 



(42) 
0. (43) 



The equation for ai(r) is longer and not really worth 
displaying. The important point is that once ao(r) is 
known, the first order equations can be solved. 

We now explore the fluctuation equations in the same 
limit assuming that the background satisfies X(t) = Xq, 
Y(t) = Yq£ 2 ^. In this paragraph we work in terms of the 
original metric variables \E', <£, x, /?, S, We only con- 
sider here the leading terms. One finds that the Newton 
potential $(r) satisfies the usual equation (recall that for 
matter c 2 = 0)) 



$ + 3H<1> + [2H + H 2 H> = 



(44) 



where 7i = a/a. (We drop the subscript o in ao because 
we work only to this order and no confusion can arise.) 
For matter a = aoT 2 and we recover, as already shown 
in section HVBl the familiar $(t) = a + b /r 5 , while for 
acceleration with a ~ l/(io — t) we have $ = ai(ro — 



All other functions are expressed via the field equations 
in terms of $ as follows: 



*(t 

(5{t 

X(r 
S(t 

/i(t 



*(t) 
2X q d 



Y§ dr 



(*(t)o(t)) 



= *(r) 
4X 



ci 



2 -^kH(r)a(r) 



-fc 2 M (r) 



C2 



Y 3 



$(t)o(t) + c 3 / dra(r) 



(45) 



where ci, 02,03 are integration constants. Recalling the 
fluid variables (|30|) it is direct to prove from here that 
both and the Newton-Poisson equation (|2l)|) are sat- 
isfied. 

The 1 1 1 series provides a systematic way to isolate the 
dark matter phase. At the same time, it provides a way 
to compute order by order deviations from ACDM which 
may reveal new features. We shall study these corrections 
elsewhere. 



V. THE COSMIC MICROWAVE BACKGROUND 
AND LARGE SCALE STRUCTURE 

The main goal in our analysis is to estimate the two 
main cosmological observables: the CMB and the large 
scale structure of the distribution of galaxies in the Uni- 
verse. 

The anisotropics in the CMB can be described in terms 
of fluctuations in temperature, (AT/T)(n) = [T(n) — 
T ] /T , where T is the average temperature in the CMB 
and T(n) is the temperature measured in the direction 
n. It is convenient to look at the variance of these fluctu- 
ations expanded in Legendre Polynomials, Pi such that 



AT AT 

;=»=-(■>'». 



E- 

1 



1 



4?r 



C/P*(n-n') 



r)+6i(r -r) 2 



where (• • • ) is the ensemble average and Ct is the angular 
power spectrum of fluctuations. We can calculate by 
evolving the Boltzman equation for the radiation distru- 
bution function, coupled to the perturbed field equations 
presented above. 

As suggested above, the evolution of EBI field as dark 
matter is exactly equivalent to that of ordinary dark mat- 
ter. Hence its effect on the CMB will be equivalent and 
we therefore expect that such observables as the peak po- 
sitions and heights will be preserved. This is clearly so if 
we look at the solid curve in Figured]- it is indistiguish- 
able from ACDM. Severe differences can arise depending 
on how the EBI field evolves in the accelerating era. If 
EBI continues to evolve as dark matter then the evolu- 
tion of the gravitational potentials is such that, again, the 
angular power spectrum is indistinguishable from that of 
ACDM. This is clearly not so in the case where the EBI 
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FIG. 4: The Cosmic Microwave Background angular power 
spectrum Ci for the models described in the realistic back- 
ground evolution. The solid curve is model A (the AEBI 
model) which is indistingushable from the best-fit WMAP-5 
ACDM model. The long-dash curve is model B, the short- 
dash curve is model C, the dot-dash curve is model D and 
the dotted curve is model E. All models have the same initial 
amplitude, same tilt (0.962) and non-zero optical depth to 
reionization (0.088) as the best-fit WMAP5 model. 

field drives acceleration. The large scale behaviour of the 
CgS is strongly dependent on the integrated Sachs Wolfe 
effect which is roughly given by: 

(^)(n) ~ £ dr(*' + tt')[(r, n(r - r)] (46) 

where primes, ', are derivatives with respect to confor- 
mal time, To is conformal time today and r* is conformal 
time at recombination. Note that the combination of po- 
tentials is the same as presented in equation (|39p and its 
evolutions is clearly different from the one experienced in 
the AEBI case. Indeed in Figure Owe plot the evolution 
of $ + W for a few cases labeled in Figure [T] Quite clearly 
the unstable, oscillatory behaviour is triggered early on 
and hence we expect it affect relatively small scales. This 
is clear from looking at Figure 2] where the modifications 
to the Cgs, through the integrated Sachs- Wolfe is present 
all the way to I ~ 150, well into the first peak. Inter- 
estingly enough, for the phantom case, the accelerating 
phase kicks in later and hence there is a smaller inte- 
grated Sachs- Wolfe effect for I > 10; once acceleration 
kicks in, however, it is much stronger than in the other 
cases and has a dramatic effect on the largest scales of 
the CfS. 

In the same way, we can directly relate the fluctuations 
in the galaxy distribution directly to the density contrast. 
It is convenient to look at the power spectrum of the 
density fluctuations by taking the Fourier transform of 



FIG. 5: The Newtonian potential combination $ + $ which is 
relevant to the integrated Sachs- Wolfe effect for the same set 
of models A-E, plotted against t/to where to is the confor- 
mal time today. In the upper panel we display model A (solid 
curve), model B (long-dash curve) and model C (short-dash 
curve). In the lower panel we show again model C (short- 
dash curve) to be compared with model D (dot-dash curve) 
and model E (dotted curve). Notice that $ + $ for models 
B-E oscillates during the transition to deSitter phase (which 
is usually a constant-w phase) , while models C-E also start to 
diverse during the deSitter phase. The presence of a bare cos- 
mological constant in model B seems to curb the divergence, 
although the oscillation remains. 

(5(k), where k is the wave number and constructing the 
variance: 

P(k) = (|5(k)| 2 ) 

Once again, we find that that the A EBI model looks 
identical to a ACDM model. In Figure [5] we plot such a 
model with a choice of parameters that render it indistin- 
guishable from the best-fit WMAP-5 ACDM model. It is 
clear that this is not true of the EBI model where we find 
that there is a strong shortage of power on small scales as 
well as a much broader turnover associated to the radia- 
tion matter transition. The effect is sufficiently dramatic 
that we don't even attempt to compare the EBI model 
to the angular power spectrum as measured by WMAP- 
5 ^28| or the power spectrum of galaxy fluctuations as 
measured by the Sloan Digital Sky Survey [29( in Figure 
[7J- the EBI model is not a viable candidate for a theory 
of structure formation- while the AEBI model is quite 
clearly a good candidate. 

The evolution of perturbations in the EBI model do 
have an interesting feature that is worth noting. In the- 
ories of pressureless dark matter, such as WIMPs, the 
evolution of perturbations is such that the two gravita- 
tional potentials are effectively identical, i.e. $ = , 3>. It 
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k (h Mpc" 1 ) k (h Mpc" 1 ) 



FIG. 6: The baryon power spectrum P(k) for the same set 
of models A-E. Once again we show model A (solid curve) 
which is indistingushable from the best-fit WMAP-5 ACDM 
model. The long-dash curve is model B, the short-dash curve 
is model C, the dot-dash curve is model D and the dotted 
curve is model E. All models have the same tilt (0.962) as the 
best-fit WMAP5 model. 



FIG. 7: The Cosmic Microwave Background angular power 
spectrum Ci (upper panel) with WMAP-5 data and baryon 
power spectrum P(k) (lower panel) for the AEBI model with 
SDSS data (model A). Both spectra are indistingushable from 
the best-fit WMAP-5 ACDM model. 



has been pointed out that in many, if not all, theories 
of modified gravity, these potentials will differ from each 
other [3Cl ] and that this may be a smoking gun for modi- 
fied theories of gravity. A plethora of observational tech- 
niques have been proposed, cross correlating galaxy sur- 
veys with weak lensing surveys and with measurements 
of the CMB [3l[ and is one of the main science targets of 
up and coming experiments such as the Euclid project. 
As pointed out in [32| such a signature is not exclusive to 
modified gravity and it suffices that the dark sector have 
a component that takes the form of anisotropic stress. 
This is indeed what we find in this theory and specifically 
in the case of EBI, where S and II can have a substantial 
effect on the evolution of perturbations. In Figure [5] we 
illustrate this fact by plotting the evolution of $ — 'J for 
a selection of models. Granted that we have been unable 
to find an EBI model that fits the data and the ques- 
tion still remains whether it is still possible to constrain 
a fundamental theory of dark matter with <f> — 'J. 



VI. DISCUSSION 

In this paper we have explored the cosmology of Uni- 
verse permeated by a field that obeys the Eddington- 
Born-Infeld equation. As shown in [9(, such a field can 
play the dual role of dark matter and dark energy and 
therefore supplies us with a counterpart to the Chaplygin 
gas as a possible unification of the dark sector. In our 
analysis we have shown that there are other regimes in 
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FIG. 8: The Newtonian potential combination $ — VP for the 
same set of models A-E, plotted against t/to where to is the 
conformal time today. In the upper panel we display model A 
(solid curve), model B (long-dash curve) and model C (short- 
dash curve) . In the lower panel we show again model C (short- 
dash curve) to be compared with model D (dot-dash curve) 
and model E (dotted curve). Notice that $ — \& for models 
B-E oscillates during the transition to deSitter phase (which 
is usually a constant-w phase) , while models C-E also start to 
diverse during the deSitter phase. The presence of a bare cos- 
mological constant in model B seems to curb the divergence, 
although the oscillation remains. 
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which the EBI field can play a different role, either as an 
alternative to simply dark matter or as a source of energy 
that can renormalize the cosmological constant. 

We have then looked at the effect the EBI field has 
on the growth of structure. We show that it can be de- 
scribed in terms of a set of fluid variables, akin to the 
construction of [27( and then identify the different key 
regimes. During the dark matter dominated regime, i.e. 
the regime in which we — and the EBI field dominates, 
the evolution of perturbation is exactly as in the stan- 
dard scenario in which the dark matter field is described 
by massive, non-relativistic particles. The gravitational 
potentials are constant and indistinguishable during this 
era, under the assumption that the initial shear and en- 
tropy is negligible. Distinctive signatures emerge in a pe- 
riod of accelerated expansion. If the EBI field dominates 
and is responsible for cosmic acceleration, there is a clear 
instability in the gravitational potentials; they not only 
grow but $ + $ diverges leading very rapidly to an over- 
whelming integrated Sachs- Wolfe effect on large scales. 
It is difficult to reconcile the angular power spectrum of 
fluctuations and the power spectrum of the galaxy dis- 
tribution predicted by an EBI theory which unifies the 
dark sector, with current data. If we restrict ourselves 
to a regime in which the EBI field simply behaves as 
dark matter, then, as expected, we find the our best fit 
model to be entirely indistinguishable from the standard, 
ACDM scenario. 

The EBI field can clearly play an important role in 
cosmology and, in particular, as a non-particulate form 
of dark matter. It's interpretation becomes interesting 
if we view the theory as bi-metric and q^ v as the true, 
geometric, metric of space-time; it is then this metric 
which is interpreted as dark matter. This is the point of 
view implicit in [8j. What we mean by the "true" metric 
of space-time is of course open to debate. Clocks and 
rulers will feel g^ v and hence real geometry will be built 
out of it and in this case q^ v plays a purely auxiliary role 
as an extra field. 

We would like to re-emphasize that the generalized EBI 
model, in which the EBI field drives cosmic acceleration, 
gives us an interesting example of a theory with an exotic 



signature: the mismatch between $ and ^ . A number 
of methods have been proposed to do ease out this mis- 
match from current and future data sets [3l|. We have 
found that the EBI field can source this mismatch with- 
out modifying gravity. This is not surprising: the EBI 
field is a two tensor with four scalar degrees of freedom. 
One linear combination of these degrees of freedom can 
be seen as anisotropic stress which can freely source the 
<f> — <!>. It turns out that its effect is severe enough that 
the integrated Sachs- Wolfe effects generated is too ex- 
treme to be reconcilable with current observations. But 
it does suggest that it may be possible to build models 
which don't modify gravity, generate accelerated expan- 
sion and could be confused with bona-fide modified theo- 
ries of gravity [32| . Consistent parameterized frameworks 
such as (33j may be able to provide alternative ways to 
distinguish such theories and it would be interesting to 
find the EBI's predictions for these frameworks. 

Finally, we would like to point out that the EBI model 
is a viable alternative to the ACDM but which may have 
particular features which make it stand out. As we saw in 
section llVBl even though the evolution of perturbations 
may be equivalent to that of that CDM if one assume no 
pressure perturbations and shear in the initial conditions, 
the non- linear evolution will be different. Pressure per- 
turbations and shear will be generated at the non-linear 
level and may play a significant role in the small scale 
structure of galaxies and clusters. Indeed, one of the ma- 
jor problems that ACDM has had to face is the excess of 
small scale power compared to observations [34| . AEBI 
may have a natural dynamical solution to this problem. 
This is one of the many aspects of this theory we wish to 
explore further. 
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